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■ Abstract. A class of linear evolutionary equations with material laws involving 

pP fractional time-derivatives is considered. The main result is well-posedness and 

causality for this problem class. The approach is illustrated with two examples: 
a fractional Fokker-Planck type equation and a class of visco-elastic materials 
described via fractional derivatives. In conclusion the possibility of imposing; initial 
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(N ! Introduction 

> 

In [6, 8j it has been shown that the standard (autonomous and linear) initial boundary value 
problems of mathematical physics share a simple common form, if considered as first order 
systems. Indeed, it is found that they are of the form 



8 V + AU = F, (1) 

where 8q denotes time-differentiation, in the usual cases A is skew-selfadjoint. The unknown 
U and V are linked by a so-called material law 

X ' 

U : V = M (d^ 1 ) U, (2) 

where Ai is a bounded operator- valued function, analytic in a ball Bq (r, r) of radius r 6 lR>o 
centered at r. A4 (Oq 1 ) is then well-defined in terms of an operator- valued function calculus 
associated with do as a normal operator in the weighted L 2 -type space H p fi (R, H), p > 
which may be described by completing the space (M, H ) of smooth //-valued functions, 
H a (complex) Hilbert space, with compact support in K with respect to the norm | • | p0 o 
induced by the inner product ( • | • )po gi ven by 



(U,V)^ [ (U(t) \V (t)) H exp(-2pt) dt. 
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We recall from [6J the Fourier-Laplace transform given as the unitary extension C p of the 
mapping 

Coo (M, H) C H pfi (R, ff) -> L 2 (R, ff) 

</? h-> / exp (—if ( • — i p)) ip (t) dt 



2tt 

is a spectral representation of the selfadjoint operator 3m do = i (<9o ~~ ip) an d so we find the 
unitary equivalence 

do = £* p (i rn Q + p)C p , 
where mo denotes the multiplication-by-argument operator given as the closure of 

Coo (R, H) C L 2 (R, ff) -> L 2 (R, iT) 
</? i y rno^p 

with 

(mo^) (A) :=Ac/?(A) in i7 

for every A G R. This observation allows us to consistently define an operator function calculus 
for operator-valued functions 



Here the linear operator Ai I im 1 , j : L 2 (R, H) — > I? (R, H) is determined uniquely by 



m +p 

for every A € R, G Coo (R, il). 

It may be worth noting that the case p G R>o utilized here leads to the forward causal integral 
representation 

{do 1 f)(t)= f f(s)ds 

J — oo 

for / G Hpfi (R, H ) as a Bochner integral. In terms of the associated operator- valued function 
calculus we also know what 

d a ,ae[o,i[, 

meana3- With this we define for 7 € R 

q^slpeq-™ (3) 



1 lt should be noted that d a is largely independent of the particular choice of p £ ]0, 00 [. Indeed, since 

r ( \ «-l 1 r ( a ) 

£ P X ]0 ,^(mo)m =—^-—^ 

we have for well-behaved ip 

P^yX ]0 ,oo[ ( m o) m o _1 * <P = do™ 1 ? 
and t 

r^jX 10 , TO[ (mo) m^- 1 * ^ (t) = ^ f ^ W *»• 
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as a natural generalization of differentiation to arbitrary real orders. The fact that (9o) 7gK is 
a family of commuting operators appears to be useful in applications, see [5] and the quoted 
literature. In contrast, the variant^ 

a Dl := dl X%>a (m) = 5™ (d^\ R>a (m)) 

with an appropriate choice of domain, known as Riemann-Liouville fractional derivative, a E R 
a parameter, |10} [3], lacks this property. This is also true for the frequently used Caputo 
fractional derivative, [TU1 13]: 

with a suitable domain. With these fractional derivatives being mere variants of limited use- 
fulness in our context, we shall utilize only the above spectral definition ([3]) for our fractional 
differentiation/integration. The aim of this paper is to inspect material laws involving frac- 
tional derivatives. For this we note that for a G [0, 1] 



?te d£ > p a (4) 

and for a G M>o 

II VII < P~ a - (5) 

Indeed, ([5]) follows from @ and to see (j3J) it suffices - by the function calculus of do - to 
consider dKt (it + p) a for t G R. We have 

<Re (it + = \it + cos (a/3) 

a/2 



p a (l + (tan/3) 2 )" cos (a/3) 
a cos (a/3) 



r (cos/3) a 

with 

/3 :=arg(it + e ) g] - tt/2, tt/2[ . 

The values of ^°^^2 are either unbounded as (3 — > ±^ (a G]0, 1[), or 1 (if a = 0, 1). In the 
first case we shall try to find a minimum. The derivative 

cos (a/3) V a (cos /3) Q sin (a/3) — cos (a/3) a (cos (/3)) Q_1 sin (/3) 

cos /3 sin (a/3) — cos (a/3) sin (/3) 



/3 i ^ a 
/3 i ^ a 



(cos/3) a+1 
sin ((a - 1) /3) 
(cos /3) 



a+l 



2 In the limit a — » — oo the spectral fractional derivative is formally recovered: 

u ~ -oo^i ■ 



3 



vanishes only for f3 = 0. Due to the behavior at the limits this must be the point where a 
minimum occurs. For f3 = we have, however, 

cos (a/3) 
(cos j3) a 

which proves (j4j. Estimates (|4]), ([5]) will be important for the discussion of the type of material 
law operators Ai (Oq 1 ) of interest here. 

Systems of the form ([1]), , give rise to equations of the form 

(d M (So" 1 ) +A)U = F, (6) 

which will be referred to as "evolutionary" equation, since the term "evolution equation" ap- 
pears to be reserved for a rather special case in this comprehensive problem class. 

The heart of theory is the positive definiteness condition: 

<Ke (U\doM (d, 1 ) U) pfifi > d {U\U) pfi>0 (7) 

holding for all sufficiently large p and some constant c\ and every U € D (do) , see e.g. [6j [7] 
for a more detailed discussion. 

In many cases of interest the material law operator has the form 

M (d l ) = Mo + d^hh + a^M 2 (d Q l ) . 

In this case condition (JJJ) is warranted if Mq is selfadjoint, for some cq G M>o 

pM + d\e Mi > c (8) 

for all sufficiently large p G M>o and 

limsup |M 2 (Oq 1 ) \ HpAR:H)eHpMRiH) < co- (9) 

The requirement ([8|) is satisfied if Mq is strictly positive definite on its range Mo [H] and 
SHeMi strictly positive definite on the null space [{0}] Mo of Mo. Here we want to consider a 
particular class of material laws of a different kind involving fractional derivatives. 

1 Material Laws with Fractional Integrals 

According to the above, requirement ([8]) is satisfied if Mo is strictly positive definite on its 
range Mo [H] and £ReMi strictly positive definite on the null space [{0}] Mo. In a little more 
detail this means that 

tSMoio : M [H] -). M [H] , k* V\t M iKq : [{0}] M [{0}] M (10) 

are strictly positive definite, where to is the canonical embedding of Mo [H] in H and kq is the 
canonical embedding of [{0}] Mo in H. Note that LqLq and kqKq are the orthogonal projectors 
onto Mo [H] and [{0}] Mo, respectively. In particular, we have 

io^o + K o^o = 1- 
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Here we want to consider material laws of a different kind involving fractional derivatives: 

M(d 1 )= ]T d^M a + d^M 2 (d^), (11) 
aenu{i} 

where 

IIC [0,1[ finite, 

M a selfadjoint and non-negative for a £ II. 

To ensure ([7]) for material law operators of the form (jlip we assume the existence of a constant 
Co € M>o such that 

^ p 1-a M a + IHe Mi > c (12) 

for all sufficiently large p E M>o and 

limsup |M 2 (0„ |^ i0(Rjfl)ffi ^ i0(KfH) < c . (13) 

Our aim is to describe condition (j 1 2[) in more explicit terms. The observation (|10p will be our 
guideline. 

With (ao, • • • , ajy) as a monotonically increasing enumeration of II we get 
p 1 ~ a M a + me Mi = ^ p 1_a *M afc + me Mi 

QGil k 

N 

= p^Ma, + Y / P 1 - ak M ctk +9teMi. 
fe=i 

We first consider the case iV = 1. Let Pq = lqLq, where to denotes the canonical embedding of 
M ao [H] in H. Then Pq is the orthogonal projector onto M ao [H] . The orthogonal projector 
onto the ortho-complement of M ao [H] , i.e. [{0}] M ao , is then 1 — Pq, which in turn is equal 
to kqKq, where kq denotes the canonical embedding of [{0}] M ao in H. We estimate 

(U\p L - a °M ao U + p 1 - ai M ai U + ^tM 1 U) = 

= (U\p 1 - a °M ao U) Q + p l - ai - {U\M ai U) Q + {U\V\tM 1 U) Q 

= p 1 -^ (t*U\ (t*M Qo t ) l* U) + p 1 -^ (L*U\L*M ai ioi* Q U) + 

+2p 1 - a ^ me (L* Q U\i* Q M ai K K* Q U) + 

V~ Ql («5^l KM ai K )K* U) Q + 

+ (t^|t^ SHeMiioiSC^o + 29 * e Wko ^cMik Kq^) + 

+ <«o^ I («S KcM iKo )k* U) . 

Absorbing the mixed terms we obtain for some eo,r/Q £ M>o 

([/l^-^M^C/ + p 1 ~^M ai U + me MiC/> = 

> ^ (iS^io) t^) + ^ (K* U\ (K* M ai K ) K* U) Q + 

+ (k* U\ (k* me MiK ) k* U) q - e (k* U\k* U) 



5 



for all sufficiently large p G R>o- Repeating this calculation we obtain with L\ as the canonical 
embedding of (kqM Qi ) [[{0}] M ao ] and K\ the canonical embedding of [{0}] (kqM Qi ko) into 
[{0}] M ao that 

(f/| / 9 1 - Q "M Q0 C/ + p 1 ~ a ^M ai U + VKzM x U) = 

> S (toM ao io) tgCOo + ^ <^£/| (^K5M ai Koti) ^£/> + 

+ {k\KqU\ (K\Kl^Kt Mikqki) k\k* U) q - e (iIk,qU\iIk%U) + 
-ei {kIkqU\k*KqU) . 

Assuming (i,Q.M ao io) j (iT'?0'^«i' s t l) an< ^ ( K i K o -/WiKo^i) strictly positive definite, we obtain 
for e\ sufficiently small as desired an estimate of the form 

(U\p 1 - a °M ao U + p 1 ~^M ai U + VKtM x U) Q > 

+ ( K >o U \ ( K *i K *o ^ M lKoKl ) k* iK * U) 

for some positive number r]\ and all sufficiently large p. The case of general JV £ N follows 
similarly by induction assuming 

L l K l-i " " " ^o M » k K o ■ ■ • Kfc-iijfc, k = 0, . . . ,N, and 
k* n ■ ■ ■ Kq 9\e M±kq • • • kjv strictly positive definite, 

where if. as the canonical embedding of • • ■ K^MaA [[{0}] M ao ] and the canonical 

embedding of [{0}] (n*^ ■ ■ ■ K$M ak K h -i ■ • • «o) into [{0}] M Qo , fc = 0, . . . , iV. 

Remark 1.1. It should be noted that the latter requirements loose much of their complexity if 
Kf, = 0, for some k < N, since then n s = for s = k, . . . , N. The simplest possible case along 
these lines is [{0}] M ao = {0} , in which M ao strictly positive definite is the only relevant 
constraint on the material law. Also note that if the above constraint (fT2|) hold for a subset 
n'cn and Mp > for j3 G n\lT then it also holds for all of IT. So, if for example M ao and Mi 
are such that (lqM^lq) and (kq^Rc M\Ko) are strictly positive definite then other coefficient 
operators Mp for /? G II \ {«o} can be left out of consideration if they are non-negative. 

Since do commutes with (<9oA4 (c^ -1 ) + A) we could extend our discussion to the Sobolev 
lattice associated with do and A + 1, see [6]. For our purposes here we do not need the 
construction of the full Sobolev lattice. It suffices to note the possibility of the ad-hoc definition 
of the dual space H„ \ (R, H) of the space H Pt i (R, H), which is the domain D (do) equipped 
with the norm | • | 10 := \8q ■ \ ■ This allows to include data such as 5® Vq, Vo G H, where 
(5 ® Vq) ($) := ($ (0) |Vb)rr for all $ G (R, H) is the appropriate representation of initial 
data in our setting. Note that these spaces give rise to the Gelfand triplet 

H p>1 (R, H) F p ,o (R, ff) ^ H Pl -i (R, H) . 

Similarly we can consider the space H Pt o (R, H\ (A + 1)), where H\ (A + 1) is the space D (A) 
equipped with the graph norm | • | 1 := \ (A + 1) • L • With H_\ (A + 1) := (A + 1) we 
obtain Gelfand triplets 

H pfi (R, H! (A + 1)) H p>0 (R, H) # p , (R, #-i (A + 1)) 

and 

Hi (A + 1) H «-»■ F_i (A + 1) . 
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Theorem 1.2. Let Ai (Oq 1 ) be a material law of the form Ull)) satisfying (TT^P. Then for every 
sufficiently large p £ R>o and every f £ H p ,k H) there is a unique solution U £ H p .k (R, H), 
k = —1,0,1, satisfying equation |6j). Moreover, the solution operator 

{d M {do 1 ) + A)' 1 : H Pyk (R, H) -> H p>k (R, if) 

is continuous and causal in the sense that, for every o£l, if g £ Li p .k (Rj H) vanishes as an 
H-valued generalized function on ] — oo, a[ then so does {8qM (d^ 1 ) + ^4) 1 g- 

Remark 1.3. 

1. Causality implies that if g G Hp,o (R, H) vanishes on ] — oo, a[ then the fractional deriva- 
tives occurring in the material law are actually Riemann-Liouville type derivatives a D^ k , 
k = 0,...,N. 

2. We have chosen to discuss the case of A skew-self adjoint in H. More general operators 
A could also be admitted by using other suitable boundary conditions, compare e.g. 
[3 HI]- The operator A can even be generalized to a monotone relation, see |llj . 



2 Some Applications 



According to Remark 11.11 the case Mq strictly positive definite is particularly simple. It has 
already been noted in [2] that such cases are easily accessible. As an application a particular 
model with Mq strictly positive definite, and IT = {0, 1/2} has been considered. Here we want 
to focus on the case Mq = 0. 

Since an additional term M2 (c^o" 1 ) satisfying (jT3j) may always be inserted into the material 
law operator (|12p . we may further simplify our perspective to the case M 2 (do" 1 ) = 0. The 
most simple case of this type is now 

II = {a} 

with a = 1 — P, P e]0, 1[, yielding 

(d^M a +M 1 +A)U = F. (15) 



2.1 Fractional Reaction-Diffusion Equations 

To be specific let us consider the fractional Fokker-Planck equation, compare e.g. [1]. This 
amounts to taking 

' V T 
V 



A :-- 



with a suitable boundary condition to ensure skew-selfadjointness, e.g. vanishing of the first 
block component $ or vanishing of the normal component of the second block component <5. 
For the material we assume 

V 1 \pio mi 
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with K a , 9\e /in selfadjoint and strictly positive definite on the respective L 2 -type block com- 
ponent spaces, /ijj. bounded, linear operators, i,k = 0,1. Thus, the system looks formally 
like 

9q a K a + fioo V T + /i Q i 

A simple row operation yields 

dl~ a n a + noo ~ (V T + noi) thi (V + /iio) 

V + filQ fin 

which reduces to a fractional convection-diffusion operator 

d^ a K a + /i 00 - (V T + noij fJ-ii (V + /Uio) = 9o _a K a + (fioo - Moi/^TiVio) - V^^V 

for the first block component, say of a solution U = I of the associated equation of 



the form (|15p . The second row represents the flux $ in terms of •& and Vi? 

$ = -/i^Vtf - MuMioi?- (16) 

Clearly for a = we recover the usual convection-diffusion equation. Speaking in terms of 
heat conduction, which is, up to a choice of units (and interpretation), governed by the same 
type of convection-diffusion equation, (fT6|) is a variant of Fourier's law. 

2.2 Fractional Visco-Elastic Material Laws 

As our second example we inspect linear visco-elasticity. It is a well-established idea that 
material laws describing visco-elastic solids may be more simple to match with measurements 
if fractional integrals are admitted, see e.g. [I], [S]. If it is assumed that M ao is strictly 
positive definite, then the lower order terms of the material law can be considered as small 
perturbations (similar to the term M2 (c^ -1 )). Thus, only an equation of the form (|15p . even 
with M\ = 0, needs to be considered in this case. The spatial operator 



A 



-Div 
Grad 



where Div is the restriction of the tensorial divergence operator div to symmetric tensors of 
order 2 and Grad is the symmetric part of the Jacobian matrix d®u of u, is skew-self adjoint if 
suitable boundary conditions on the boundary of the underlying domain £1 C M 3 are imposed, 
e.g. vanishing of the displacement u on the boundary. This choice of boundary condition 
would amount to replacing Grad by the closure Grad of the restriction of Grad to vector 
fields with smooth components vanishing outside of a compact subset of Q. The above general 
considerations allow for very general materials to be covered by the approach. 

Let us take a closer look at a fractional Kelvin- Voigt model for visco-elastic material. In this 
model we have the equation 

r]dlu - DivT = / 



•s 



linking the stress tensor field T with the displacement vector field u, accompanied by a material 
relation of the form 

T=(C + Dd$)£ (17) 

with £ := Grad u, a £ [0, 1] . Here the case D = would correspond to purely elastic behavior 
and the case a = 1 and D ^ would describe the classical Kelvin- Voigt material. Introducing 
v := 8qu as a new unknown we arrive - by differentiating (|17|) - at 

7]dov — DivT = / 
d {C + Dd$y 1 T = Gradu 

yielding formally an evolutionary equation of the form 

Po \ / - Div \ \ / v\ ( f 



do \0(C + Dd$)- 1 ) + y-Grad J J \T J \0 

If C, D both are assumed to be selfadjoint and D is strictly positive we have for a E]^, 1] in 
the above terminology 



with 



Po 




and 



M 2 := X {do 1 ) - M - M a d^ a 



Remark 2.1. A more careful inspection shows that the material can be more general. It suffices 
to assume a E]i, 1], C,D > and 

C + Dp a >c (18) 

for some cq E M>o and all sufficiently large p E M>o- Estimate (fT8|) holds if 7Tq D Ciro t £>, 
ir^ D Diri t D are both strictly positive definite, where tvq d an d ki d are t ne canonical embed- 
dings of null space and range of D in the underlying Hilbert space of symmetric 3 x 3-matrices 
with entries in L 2 (Q), respectively. 



2.3 Initial Boundary Value Problems 

The question may be raised, how initial value problems could be posed for problems of the 
form (|15p . Since fractional derivatives of order f3 E]0, 1[ are non-local, i.e. we have memory 
effects, prescribing a pre-history appears to be more appropriate. This amounts to given 
Hp$ (R, iT)-data prior to the initial time 0, i.e. "the pre-history", and is clearly covered by the 
above. 

However, if one would force the issue of initial conditions, there appear to be two different ideas 
to implement initial conditions. The first option would be to consider with F E Ft p ,o QSL,H), 
W a E M a [H] , 

Q^MaU + M X U + AU = F + 5 (8) W a (19) 
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as a straight-forward generalization of the case /3 = 1, see [6]. This translates to 

d do a M a U + M\U + AU = F + 5®W a , 

which following the rationale in [6] would be in point-wise terms an implementation of the 
initial jump condition 

[d a M a U - XtL>0 ® W a ) (0+) - (dv a M a U - Xr>0 ® W a ) (0-) = in 

If F = on R<o then causality yields 

(d^MJJ - Xr>0 ® W a ) (0-) = in 

and so, in this case an initial condition of the form 

(d^ a M a U) (0+) = W a in fr_! 

results. 

Alternatively, let 

and G 6 H p ^ (R, i7) be given and consider 

0% [M a V - Xs>0 ® V a ) + MtV + AV = G. (20) 

For the case a = 0, f3 = 1, G = F and Vo = Wo this would be a re- formulation of (|19p . This 
justifies to take (|20|) as an appropriate alternative generalization. Problem (I20|) is equivalent 
to the evolutionary problem 

d?M a V + MxV +AV=G+ d$Xu >0 ® ^ 

= G + a - Q 5®y a . 

Both approaches are clearly closely related. Indeed, for G = dQ a F and V a = W a we see by 
comparison that 

v = d^ a u. 
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